Magnetization plateaus for spin-one bosons in optical lattices: Stern-Gerlach 
experiments with strongly correlated atoms 
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We consider insulating states of spin-one bosons in optical lattices in the presence of a weak mag- 
netic field. For the states with more than one atom per lattice site we find a series of quantum phase 
transitions between states with fixed magnetization and a canted nematic phase. In the presence of 
a global confining potential, this unusual phase diagram leads to several novel phenomena, including 
formation of magnetization plateaus. We discuss how these effects can be observed using spatially 
resolved density measurements. 



Far-off-resonance optical traps can confine neutral 
atoms regardless of their hyperfine state and open 
exciting possibilities for exploring multicomponent sys- 
tems of degenerate atoms |ll!SI!E0HSHIll 
El EE E1E1 E3- Experimental studies of the hyperfine 
spin-one manifold of 23 Na atoms confined in a single op- 
tical trap already revealed the grou nd state spin struc- 
ture in external magnetic field |l8j . the formation and 
persistence of the metastable spin domain configurations 
[19j , and quantum tunnelling across spin domain bound- 
aries [20| . These phenomena have been understood using 
mean- field theory [HHEl since in a large trap atoms can 
be considered as weakly interacting. At the same time, 
theoretically predicted many-body features of spin-one 
condensates, such as a spin singlet nature of the ground 
state for even number of atoms |4j , could not be observed 
in such experiments, since the energy difference between 
various spin states is inversely proportional to the volume 
of the system and is extremely small for realistic traps 



Several approaches have recently been suggested to 
creating strongly correlated quantum states of spinor 
atoms in optical lattices produced by standing wave laser 
fields [3, la la LLii LLJ| • In the presence of a deep optical lat- 
tice wavefunctions of atoms are localized near the lattice 
potential minima, which leads to a strong enhancement 
of interactions between atoms |22t I23L l24j and can result 
in dramatic changes in the properties of ultracold gases. 
It is also important that in optical lattices one can have 
a small number of atoms per lattice site (in experiments 
of Ref. |23| this number was around 1-3). In an insulat- 
ing state hopping between different sites is negligible and 
each well behaves essentially as a small condensate. In 
this regime, the behavior of spin-one bosons in each well 
will be dominated by the strong spin-spin interactions 
such as those predicted by 4] . Two component atom 
mixtures in optical lattices have already been realized in 
experiments of Ref. [2(| . 

Non-trivial nature of the ground state of localized spin- 
one bosons in optical lattices has been discussed previ- 
ously in 0, 0, E| • Virtual hoppings of the atoms give rise 
to spin exchange interactions between neighboring sites 



and lead to several distinct insulating phases, which differ 
in their spin correlations. In two and three dimensional 
lattices states with odd numbers of particles per site are 
always nematic, and for states with even numbers of par- 
ticles per site, there is always a spin singlet phase, and 
there may also be a first order transition into the nematic 
phase. In this paper we extend earlier analysis and show 
that in a magnetic field insulating states with more than 
one atom per site undergo a series of phase transitions 
between spin gapped phases with quantized magnetiza- 
tion and phases with so called canted nematic order, in 
which magnetization can vary continuously. The critical 
properties of these phase transitions have been recently 
investigated in [15J . Magnetic phase diagram is similar to 
the phase diagram of the spinless boson Hubbard model: 
states with quantized magnetization are analogous to the 
Mott insulating phases, while the canted nematic phase 
is analogous to the superfluid phase. We also propose two 
kinds of experiments that can verify the presence of such 
magnetization plateaus and demonstrate the many-body 
nature of insulating states of spin-one bosons in optical 
lattices. 

An effective Hamiltonian for spin-one bosons in an op- 
tical lattice in the presence of a magnetic field is given 



by i 



K = -t (4r°j* + a] a a ll7 ) + -J- ^ fii[hi - 1) 



i i i 

Here a\ a are creation operators for particles in the low- 
est Bloch band localized on site i and having spin 
components a = 0, ±1; rij = a\ a ai a and Si — 
J2aa' ^ia^oa'O-ia 1 are the number and spin operators for 
site i {T aa i are the usual spin matrices for spin 1 parti- 
cles). For each well the collective spin of the atoms satis- 
fies constraints Si < Ni and Si + Ni is even. Parameters 
t, Uq and U2 for a realistic case of a three dimensional 
cubic lattice have been obtained in ||. The ratio of the 
interaction terms in Q, U2/U0, is fixed by the ratio of 
the scattering lengths and is independent of the nature of 
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the lattice. Scattering lengths of 23 Na obtained in Ref. 
[2^ give U2/U0 = 0.04. In this paper we neglect effects of 
the quadratic Zeeman shift since magnetization plateaus 
that we are interested in appear for magnetic fields of the 
order of mGauss (assuming a typical J7 2 ~ 0.1kHz). 

When the spin dependent interaction (U2) is much 
smaller than the Hubbard repulsion (Uq), the superfluid 
- insulator transition |24],|2j| is determined mostly by Uq. 
The spin gap U2 term, however, is important inside the 
insulating regime, in which nontrivial spin phases appear 
as a result of a competition between a spin gap, a mag- 
netic field, and spin exchange interactions, induced by 
fluctuations in the particle number. The magnetic phase 
diagram can be most easily understood by considering 
the limit of large number of atoms per site, N ^> 1. In 
this case the effective spin Hamiltonian for the insulating 
state can be written as a model of quantum rotors, inter- 
acting via rotationally invariant quadrupolar interaction 
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Angular momentum Si 
lective spin on site i and inherits constraints of the micro- 
scopic Hamiltonian . In the mean field approximation 
[30L l3Ij we replace by a sum of single site Hamiltoni- 
ans 



H 
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zN 2 t 2 1 
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with z being the number of nearest neighbors, and impose 
a self-consistency condition on the nematic order param- 



eter Q ab = {n ia nib} 



Magnetic phase diagrams 



obtained by solving self-consistently for the cases of 
odd and even filling factors are shown in Figs. ^and[21 (for 
details see When t = the Hamiltonian © does 

not mix different spin eigenstates. In a magnetic field the 
system has a series of level crossings between states with 
different values of the spin. Each of these states has a 
gap in the excitation spectrum and remains stable after 
turning on a finite value of t. This results in lobes of fixed 
magnetization shown in Figs. ^ and [21 Only at points 
where spin eigenstates come into degeneracy exchange 
interactions give rise to mixing of different spin eigen- 
states, leading to a Canted Nematic phase. The latter 
has an expectation value of the nematic order parameter 
Qab in a plane perpendicular to the magnetic field 
Hence, it has spontaneous breaking of the symmetry of 
spin rotations around the direction of magnetic field. For 
sufficiently large magnetic field the system becomes fully 
polarized with Si = N. We note that for N = 1 and 
N = 2 one can derive effective spin interactions that do 
not rely on the large N approximation of equation (0) 
[j| . These models give magnetic phase diagrams that are 
qualitatively similar to Figs. n an d|21 28]. We also point 




FIG. 1: Magnetic insulating phase of spin-one bosons in opti- 
cal lattice with an odd number of atoms per site N. Magneti- 
zation is fixed inside the lobes but varies continuously inside 
the Canted Nematic phase. The latter has an expectation 
value of the nematic order parameter Q a ^ in the plane per- 
pendicular to the direction of magnetic field. For sufficiently 
large magnetic field the system becomes fully polarized with 
S = AT (not shown here). The insert shows the ground state 
magnetization (per site) as a function of magnetic field for 
2zN 2 t 2 /Uq = 5C/2- Note, that near H = magnetization 
grows linearly with magnetic field and quickly reaches the 
Ssite.z = 1 plateau. 
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FIG. 2: Magnetic insulating phase in optical lattice with an 
even number of atoms per site. The insert is for 2zN 2 t 2 /Uo = 
2U 2 . 



out that inside the lobes of fixed magnetization and in the 
mean-field approximation the many-body wavefunctions 
factorize |¥) = Hi \Si = S, S lz = S). 

A typical experimental system has a parabolic con- 
fining potential in addition to the periodic optical lat- 
tice. This means that the chemical potential is not uni- 
form and for small t the system breaks into insulating 
domains with different integer filling factors, separated 
by regions of the superfluid phase (see Refs. 0> 
for analysis of the spinless case) . Another important ex- 
perimental constraint is that the total magnetization is 
fixed by the initial state of the system. In the discus- 
sion above we showed, that when we change magnetic 



field, it is energetically favorable to adjust the magneti- 
zation in order to utilize some of the Zeeman energy (see 
Figs, ^and (2J . From equation JTJ we observe, however, 
that the spin component parallel to the applied field is 
conserved. For example, if the magnetic field is along 
the z axis, tSt t,z — ^iz is ^ good quantum number 
of the system (spin non-conserving interactions, such as 
the dipolar relaxation, are typically small) and should 
not change even when H z is changing. For a single large 
trap this feature allowed to study magnetic properties 
of spinor condensate 01 a t magnetic fields for which the 
true ground state should be fully polarized. 

A way around spin conservation has been demon- 
strated in and relies on applying spatially varying 
magnetic fields and performing Stern-Gerlach imaging. 
Here we extend these ideas and suggest an approach to 
experimental observation of the magnetization plateaus 
discussed above. The idea of our first experiment is 
shown in Fig. [3] We consider a strongly anisotropic 
trap in which a magnetic field gradient is applied paral- 
lel to the long axis of the system. When there is a mag- 
netic field gradient parallel to the long axis, there should 
also be gradients in the transverse directions (both VH 
and V x H should be zero). We assume that the size 
of the condensate in the transverse directions is small 
enough that we can neglect the effects of the magnetic 
field in transverse directions. However, the size of the 
condensate should be larger than the optical lattice pe- 
riod in any direction, so that we can consider it as a 
three dimensional system. As we discussed before, the 
uniform part of the magnetic field has no effect on the 
state of the system, so in our discussion we set it to 
zero. To be concrete, we assume that the largest in- 
sulating domain at the center of the trap has a filling 
factor N — 6 and that our system has been prepared 
to have Stot.z = 0. In a nonuniform magnetic field dif- 
ferent parts of the N = 6 domain minimize their en- 
ergy for different values of magnetization. When the 
field gradient is sufficiently large and tunnelling is small 
enough, the locally favored magnetization changes in a 
step like fashion from 5 s it e ,z = — 6 (per site) on the left 
to iSsit^z = 6 on the right. Such state is also consistent 
with the spin conservation, since St t,z remains zero. So, 
the configuration that minimizes the energy has plateaus 
in the spatial profile of magnetization, with spin polariza- 
tions (per atom) Satom.z = -1,-2/3,-1/3,0,1/3,2/3,1, 
as shown in Fig. |31 Each plateau has a length of the 
order of 2£/2/|V-ff|. The appearance of spin plateaus 
in a nonuniform magnetic field is analogous to the do- 
main structure of condensates in optical lattices in the 
presence of a non-uniform global confining potential that 
was discussed before 0, . In the latter case the to- 
tal density is fixed by the number of atoms in the trap, 
but insulating phases with different integer filling fac- 
tors exist due to the confining potential. To detect the 
magnetization plateaus shown in Fig. [21 one needs to 
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FIG. 3: A system of spin-one bosons in an optical lattice 
confined by a parabolic potential. In the insulating regime 
the cloud breaks into insulating domains with different inte- 
ger filling factors. When magnetic field gradient is applied 
parallel to the long axis of the trap, magnetization plateaus 
develop inside individual insulating domains. 



image different parts of the trap separately, measuring 
spin polarization per atom as a function of the position: 
(S'atom.z) = (n+ - n-)/(n + + n + U-), where n±,n 
are the local densities of atoms with a = ±1 and re- 
spectively. This quantity can be most easily measured 
in Stern-Gerlach time-of-flight experiments. If a small 
gradient of magnetic field is applied during expansion, 
clouds with different spin components spatially separate, 
and one can measure the number of atoms with different 
spin components using light scattering. 

To illustrate system parameters needed to realize this 
experiment we consider a cigar shaped condensate of sizes 
400 x 10 x 10 jjm. For an optical lattice created with 
A = 985 nm lasers and 4 x 10 5 atoms in a condensate 
we get the maximum density of six atoms per well in the 
center of the trap. To observe five plateaus in this setup, 
one would need magnetic field gradients ~ lOOmG/cm, 
and spatial resolution of 30/im. These parameters have 
already been achieved in experiments of Ref. [18(. In- 
side the insulating phase the characteristic timescale for 
spin relaxation between different wells is set by the ex- 
change interactions r ex = TiUq/ (Nt) 2 . For t = O.lkhz 
and Uq = 2khz we find times of the order of hundreds of 
milliseconds. So, if a magnetic field gradient is applied in 
the insulating regime, one needs to wait at least that long 
for magnetic plateaus to develop. Experimentally it may 
be more efficient to apply magnetic field gradient when 
the system is in the superfluid regime and then take the 
system to the insulating state by slowly reducing t. 

We now discuss our second " spin decoration" approach 
to experimental detection of spin gap effects in the insu- 
lating regime of spin-one bosons. The idea of this method 
is that when a system is prepared with a non-zero magne- 
tization in the absence of magnetic field gradients, mag- 
netization gets distributed non-uniformly among insulat- 
ing domains with different filling factors. Spin polariza- 
tion appears predominantly in domains with odd num- 
bers of atoms per site and is pushed out of the regions 
with even fillings. To justify this conjecture we propose 
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FIG. 4: A system of spin-one bosons with a nonzero to- 
tal magnetization in an optical lattice (no external magnetic 
field). Magnetization gets distributed among regions with 
odd filling factors and is pushed out of the regions with even 
fillings. 



the following argument. In the regime, when tunneling 
is smal\(zN 2 t 2 /Uq <C U2), there is a crucial difference 
in spin susceptibility between odd and even phases. For 
even domains, one needs to pay an energy SU2 to break a 
singlet state and have <S s itc,z ~ 2. For odd sites the lowest 
energy state already has 5 = 1, and energy cost to polar- 
ize existing spins is of the order of zN 2 t 2 /Uq. Therefore, 
externally imposed nonzero magnetization will be redis- 
tributed in odd insulating domains. For small magneti- 
zation per site, energy goes as 



1 



2Xodd(A0 



S' 2 



(4) 



If Xodd(^V) was the same for all N, then magnetization 
would be distributed uniformly among all odd domains. 
In reality, Xodd(N) decreases with increasing N, so we 
expect larger magnetization for insulating domains with 
smaller number of atoms. Quadratic dependence in l|4*|l 
ensures, however, that all domains with odd filing fac- 
tors acquire finite magnetization. So, in experiments we 
expect to find a picture of alternating even and odd do- 
mains, in which odd domains have finite magnetization 
and even domains have none (see Fig. This picture is 
valid until all odd regions have magnetization S'sitcz = 1- 
For the experimental setup discussed earlier this corre- 
sponds to (S a tom,z) i; 0.1. Performing spatially resolved 
measurements of spin polarization one should be able to 
observe such a modulated structure of magnetization. 

In summary, we discussed insulating states of spin-one 
bosons in optical lattices in the presence of a magnetic 
field. For systems with more than one particle per site we 
demonstrated the existence of a series of phase transitions 
between phases with fixed magnetization and the canted 



nematic phase in which magnetization can vary contin- 
uously. We considered experimental signatures of this 
novel magnetic phase diagram, including formation of 
magnetization plateaus in the presence of a magnetic field 
gradient and a modulated spin density in nonuniform sys- 
tems with non-zero total magnetization. We thank E. 
Altaian, D. Petrov, A. Sorensen, DW. Wang, and A. 
Polkovnikov for useful discussions. This work was par- 
tially supported by NSF (PHY-0134776,DMR-0132874), 
Sloan and Packard Foundations. 
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